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ABSTRACT 

We present a model for gamma-ray bursts (GRB's) in which a stellar mass 
black hole acquires a massive accretion disk by merging with the helium core of 
its red giant companion. The black hole enters the helium core after it, or its 
neutron star progenitor, first experiences a common envelope phase that carries 
it inwards through the hydrogen envelope. Accretion of the last several solar 
masses of helium occurs on a time scale of roughly a minute and provides a 
neutrino luminosity of approximately 10 51 - 10 52 erg s _1 . Neutrino annihilation, 
0.01% to 0.1% efficient, along the rotational axis then gives a baryon loaded 
fireball of electron-positron pairs and radiation (about 10 50 erg total) whose 
beaming and relativistic interaction with circumstellar material makes the GRB 
(e.g., Rees & Meszaros 1992). The useful energy can be greatly increased if 
energy can be extracted from the rotational energy of the black hole by magnetic 
interaction with the disk. Such events should occur at a rate comparable to 
that of merging neutron stars and black hole neutron star pairs and may be 
responsible for long complex GRB's, but not short hard ones. 



Subject headings: gamma rays: bursts, gamma rays: theory, black hole physics 



1. INTRODUCTION 

Many close, massive binaries are expected to pass through an evoutionary phase during 
which a compact object (neutron star or black hole) enters the envelope of its red giant 
companion. Some of these systems eject the hydrogen envelope and eventually evolve into 
doubly degenerate binaries such as the Hulse- Taylor pulsar system (Hulse & Taylor 1975, 
van den Heuvel 1995). In others, however, the orbital energy of the system is insufficient 
to remove the hydrogen envelope and thus shut off drag before the compact object enters 
the helium core. As the compact object and helium core coalesce, the helium core is tidally 



- 2 - 



disrupted into an accretion disk around the compact object with a radius equal to a fraction 
of the initial helium core, ~10 9 — 10 10 cm. The accretion rate onto the compact object 
may be initially as large as the Bondi-Hoyle rate, almost a solar mass per second (Houck & 
Chevalier 1991; Brown 1995; Chevalier 1996; Fryer, Benz, & Herant 1996). Later, following 
disk formation, the accretion rate is limited by the viscous time scale, (afl)' 1 ~ 100 s, 
where a is the disk viscosity parameter, typically 0.1, and Q is the orbital angular velocity. 
During the merger process, the compact object becomes a black hole, if it was not already. 

The resultant object, a ~ 3 M Q black hole accreting the remaining several solar masses 
of helium, resembles, except for the larger store of available mass and higher specific 
angular momentum, j ~ 10 18 cm 2 s _1 , the class of accretion-disk/black- hole GRB models 
based on the merger of neutron star pairs (NS/NS) or black hole neutron star (BH/NS) 
pairs. It also resembles, except for its larger angular momentum, rotating failed supernovae 
(Woosley 1993). Because of the large angular momentum, the accreting black hole will 
likely become a Kerr hole with a rotational speed near the maximum allowed by relativity. 
As the black hole spins up, the efficiency at which the gravitational energy is released 
rises from 0.06 for Schwarzschild (non-rotating) geometries to as much as 0.42 Mc 2 for 
Kerr geometries. Taking 0.1 as a representative value of this efficiency, the accreting black 
hole's neutrino luminosity for 100 seconds is about 10 51 — 10 52 erg s _1 . Conversion of 0.01 
to 0.1% of this energy (see below) into a pair fireball by neutrino annihilation along the 
rotational axis leads to burst energies in the range of 10 49 - 10 51 erg, equal to those achieved 
by NS/NS mergers (Ruffert et al. 1997). The strong beaming that will likely accompany 
accretion-disk/black-hole GRB models then allows sufficiently-high effective energies to 
power the GRB. 

As we shall also see, the rate of such helium mergers may be an order of magnitude 
greater than that of NS/NS and BH/NS binaries combined. GRB's from helium mergers 
differ from those of NS/NS and BH/NS mergers in that no appreciable gravitational wave 
signal should precede a helium merger gamma-ray burst as they would for gamma-ray 
bursts from the merger of NS/NS or BH/NS binaries. In addition, gamma-ray bursts from 
helium mergers occur immediately at their star-formation sites, whereas the delay in the 
merger of NS/NS and BH/NS binaries allows these gamma-ray burst progenitors to leave 
not only the site of their formation before merging, but, in some cases, their host galaxy as 
well. 
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2. 



PRE-BURST EVOLUTION 



Our model begins with a close binary system of two massive stars (each greater than 
8 M Q ). As the more massive star (primary) evolves off the main sequence, it fills its 
Roche- lobe and transfers mass onto its companion (secondary). If the system remains 
bound after the supernova explosion of the primary, a binary composed of a neutron 
star and a massive main-sequence star results, possibly a massive X-ray binary. When 
the secondary, in turn, evolves off the main sequence and overflows its Roche lobe, the 
neutron star enters a common envelope with the secondary and begins to spiral towards the 
secondary's helium core. The inward motion continues for roughly 1 (Sandquist et al. 1998) 
to 1000 years (Taam, Bodenheimer, & Ostriker 1978) until either the secondary's hydrogen 
envelope is completely ejected and the system becomes a helium star/neutron star binary 
or the compact star merges with the helium core. 

Whether or not the neutron star then finally merges with its secondary companion 
depends upon one of the most uncertain parameters in population synthesis studies: the 
common envelope efficiency (qce)- This single parameter describes the efficiency at which 
the orbital energy, when injected into the companion's hydrogen atmosphere as the neutron 
star inspirals during a common envelope phase, drives off the companion's envelope. Using 
o;ce, on e can estimate the final orbital separation (Af) in terms of the initial separation 
(Ai) of the binary after common envelope evolution (Webbink 1984): 



where M s , M# e , and M^s are, respectively, the masses of the secondary, the secondary's 
helium core, and of the neutron star and = Ri/Ai is the dimensionless Roche lobe radius 
of the secondary. Some hydrodynamical models have simulated this inspiral, but only 
for specific systems. Although a general consensus of the value of the common envelope 
efficiency has not been achieved, best estimates give c*ce ~ 0.5 (Taam et al. 1997). 

During the inspiral, the neutron star accretes material at the Bondi-Hoyle rate, 
releasing the accretion energy via neutrino emission (Houck & Chevalier 1991, Chevalier 
1993, Chevalier 1996, Brown 1996; Fryer, Benz, & Herant 1996). Bethe & Brown (1998) 
have calculated that in the hydrogen inspiral alone, a neutron star is likeley to accrete 
~ 1M , leading to the collapse of that neutron star into a black hole. In the tenuous layers 
of hydrogen envelope, Chevalier (1996) showed that if the infalling material had sufficient 
angular momentum, the temperature of the material on the neutron star surface may not get 
high enough to emit neutrinos, and the Eddington limit would constrain the accretion rate. 
However, near or below the surface of the helium star, where most of the accretion occurs, 
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the densities (and the subsequent Bondi-Hoyle accretion rates) are so high > 1M /yr, 
that angular momentum will not prevent super-Eddington accretion via neutrino emission 
(Chevalier 1996). The analysis of Chevalier assume that all of the mass accretes onto 
the neutron star flows through the disk, but accretion along the angular momentum axis 
dominates the flow (Fryer & Kalogera 1998). In addition, 3-dimensional simulations of 
Bondi-Hoyle accretion suggest that the amount of angular momentum accreted accreted 
in common envelopes may be much less than that assumed by Chevalier (e.g. Ruffert & 
Anzer 1995, Ruffert 1997), again weakening the effects of angular momentum. Thus, for all 
of these reasons, the Bethe & Brown (1998) estimate of the accreted mass is valid. 

As the compact object merges with the helium core, the Bondi-Hoyle accretion rate in 
the helium core can reach ~ lM Q /s. From population synthesis calculations using the code 
developed by Fryer, Burrows, & Benz (1998), the helium core and compact object have 
average masses of ~ 4M ,2M respectively at the time of merger, and the mean specific 



angular momentum of the system is j m MbhMhcJ GRn c / (M B}1 + M Hc ) 3 ~ 10 18 cm 2 /s 



where G is the gravitational constant and Mbh is the black hole mass when it reaches the 
helium core radius (Rn c ). As the two objects coalesce, the orbital energy will drive off 
what remains of the hydrogen envelope as well as some of the helium core. The angular 
momentum will be injected into the system to form a rotating disk of helium around the 
black hole. 



Whether or not the resultant system is a feasible model for gamma-ray bursts depends 
upon the accretion rate of the disk onto the black hole. Based upon the vv annihilation 
paradigm for black hole-disk GRB models (Goodman, Dar, & Nussinov 1987; Paczynski 
1991;; Woosley 1993; Ruffert et al. 1997), we can estimate the energy produced in 
electron-positron pairs for our model. For thin disk accretion, the efficiency at which 
gravitational potential energy is emitted in radiation (neutrinos) is well known (e.g, Shapiro 
& Teukolsky 1983): 5.7% for Schwarzchild black holes, and 42.3% for maximally rotating 
Kerr black holes. The black hole not only gains mass from the accreting material, but also 
angular momentum, roughly equal to the angular momentum of the material at the last 
stable orbit. The black hole will reach its maximal rotation after accreting AM = 1.846Mbh 
for thin-disk accretion (Thorne 1974), or AM « Mbh for thick-disk accretion (Abramowicz 
& Lasota 1980). In the helium merger model, ~ 4 — 5M accrete onto the 1.4M progenitor 
neutron star/black hole. The early accretion will accrete along the angular momentum 
axis, but enough may accrete along the equator to spin-up the black hole, and gravitational 
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energy conversion efficiencies may be as high as 10-20% for these models. 

This neutrino emission, if emitted in a disk, is then converted into electron-positron 
pairs with an efficiency (Ruffert et al. 1997): 

where L u is the total neutrino luminosity (all flavors), < e v > is the mean neutrino energy, 
and Rd is the inner disk radius. To estimate a maximum energy in electron-positron pairs 
which will then drive the GRB, we assume an accretion rate equal to the Bondi-Hoyle rate 
(lM /s) for a Kerr black- hole accreting 1M Q : ~ 10 52 erg. However, for a Schwarzchild 
black hole accreting 1M at accretion disk rates of M^^/ (aVt) = O.OlM /s, the energy 
drops to 10 48 erg. Typical values for a 4M accretion disk scenario lead to energies between 
10 49 — 10 50 erg, quite comparable to the burst energies of NS/NS mergers (Ruffert et al. 
1997). 

A more promising model relies upon a strong magnetic field being produced in the 
accretion disk, which can then tap the rotational energy of the black hole to power a GRB 
(Blandford & Znajek 1977; MacDonald et al. 1986; Paczynski 1991, 1997; Woosley 1993; 
Rees & Meszaros 1997, Katz 1997): 

. / jc \ 2 /M BH \ 2 / B \ 2 _ , 



where j is the specific angular momentum of the black hole and B is the magnetic field 
strength in the disk. The large disks of our helium merger model will spin-up the black 
hole, providing sufficient energy and power for strongly magnetic fields to produce more 
than 10 52 erg of beamed energy over 100 s provided that a fraction of order 1% of the 
equipartition field value, B 2 /8n ~ pv 2 is attained in the inner disk. 



4. EVENT RATES 

Our burst energy estimates are well within limits to drive a GRB, especially if the 
strong beaming, which is likely to occur, is included. But for helium mergers to be a 
viable gamma-ray burst model, they must also have a sufficiently high formation rate to 
explain the observations. Population synthesis studies of massive binaries is fraught with a 
variety of unknown parameters: e.g., kicks imparted onto neutron stars at birth, common 
envelope efficiency (ckce), the initial mass function, the binary mass ratio distribution. 
Even the stellar radii during giant phases are not known to accuracies better than a factor 
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of ~ 2 — 5. Using a slightly modified version of the Monte Carlo code described in Fryer, 
Burrows, & Benz (1998), we have run a series of population synthesis calculations (Fryer, 
Woosley, & Hartmann 1998). Here we present specific results of two simulations using two 
delta-function kick magnitudes (50,150 km s" 1 ) directed isotropicallyQ. 

In Figure 1, we compare population synthesis results of the formation rate vs. age 
of gamma-ray bursts from helium-mergers with that from NS/NS and BH/NS binaries 
combined for a galaxy with a burst of star formation and a galaxy with a constant supernova 
rate of 10 -2 y -1 . These rates depend sensitively on the supernova kick and many of the 
binary parameters which may alter the formation rate by over an order of magnitude (Fryer, 
Burrows, & Benz 1998; Fryer, Woosley, & Hartmann 1998), but the rate for helium mergers 
remains comparable to, and generally greater (often by an order of magnitude) than, the 
gamma-ray burst rate from NS/NS and BH/NS binaries combined. The rate is enough to 
provide the observed bursts even if a large beaming factor is invoked (Wijers et al. 1998). 

Whereas helium merger GRB's tend to occur in the star formation regions in which 
they are born, NS/NS and BH/NS binaries may not merge until long after their creation 
and may, therefore, leave these regions before exploding as a GRB. Many NS/NS and 
BH/NS binaries form with systemic velocities which not only drive them out of their places 
of birth, but may also drive them beyond their host galaxies. Assuming that the host 
galaxy's gravitational potential has no effect on the outmoving binary (systemic velocities 
are ~ 100 — 200 km s -1 , so this assumption is valid for low-mass galaxies), we can estimate 
the distribution of distances these binaries travel prior to merging as a GRB (Figure 2). 

The dominant factor governing the distance estimates is the merger timescale for the 
binaries. This timescale, in turn, depends sensitively upon the relation of the helium star 
radius with its mass. For our simulations we use following helium star mass-radius relation^: 

f 2.398 - 2.013 log M He M He < 2.5M 

g He ' max \ -0.699 + 0.0557 (log Mhc-0.172)- 2 - 5 M He > 2.5M . [ ) 

For those compact objects whose inspiral take them within the helium star radius, the 
system merges and becomes a helium merger GRB, not a NS/NS or BH/NS binary. Hence, 
to form a NS/NS or BH/NS binary, before the supernova explosion of the helium star, 
the helium star/compact object separation must be greater than this radius. This limiting 



1 We set cuce = 0.5, we use a Scalo (1986) initial mass function and a flat mass ratio distribution. See 
Fryer, Burrows, & Benz (1998) and Fryer, Woosley, & Hartmann (1998) for details. 

2 This relation uses a radius scaled down by a factor of 5 from the mass-radius relation derived in Kalogera 
& Webbink (1998) to match the radii from Woosley, Langer, & Weaver (1995) 
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pre-explosion separation defines the orbital separation distribution after the supernova 
explosion of the helium star. As the helium star radius decreases, so then does the mean 
orbital separations of the NS/NS and BH/NS binaries, and hence, the merger times and 
distances traveled before merger decrease as well. In figure 2, we see that the mean distance 
traveled at the time of merger is roughly ~ 10 — 100 kpc. If the helium star radii were a 
factor of 5 higher, this mean would increase to ~ 1 Mpc. Thus, it is possible that GRB's 
driven by the merger of compact binaries may not occur in their host galaxy, whereas 
helium mergers will all occur therein. 

The authors appreciate stimulating conversations with Hans Bethe and Gerry Brown 
regarding common envelope evolution and black hole formation. We'd also like to thank 
Peter Bodenheimer for his advice on neutron star inspiral calculations, Bob Popham for 
sharing with us output from his accretion disk calculations, and Alex Heger for helpful 
suggestions. This work was supported by NASA (NAG5 2843 and MIT SC A292701) and 
the NSF (AST 94-17161). 
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Fig. 1. — (a) The number of gamma-ray bursts vs. age for a single burst of star formation 
with 10 8 supernovae. Gamma-ray bursts continue to occur from mergers of NS/NS and 
BH/NS binaries long after the initial burst of star formation. Helium merger gamma- 
ray bursts occur along with the star formation and should not be observed in old stellar 
populations, (b) The number of gamma-ray bursts produced per Myr vs. age assuming 
a constant supernova rate of 0.01 y _1 . The helium merger gamma-ray burst rate quickly 
reaches a peak and then remains flat. 
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Fig. 2. — Gamma-Ray Burst Location. Distribution of gamma-ray bursts vs. distance from 
formation site assuming no effects from the galactic potential for the NS/NS and BH/NS 
binary mechanism. Helium merger gamma-ray bursts will explode in the star formation 
region in which they were produced whereas some gamma-ray bursts from degenerate binaries 
may extend as far as 1 Mpc from their host galaxy. 



